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Steady-state visual evoked potentials and phase synchronization in migraine 
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We investigate phase synchronization in EEG recordings from migraine patients. We use the 
analytic signal technique, based on the Hilbert transform, and find that migraine brains are char- 
acterized by enhanced alpha band phase synchronization in presence of visual stimuli. Our findings 
show that migraine patients have an overactive regulatory mechanism that renders them more sen- 
sitive to external stimuli. 
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Phase synchronization was introduced for coupled 
chaotic systems by Rosenblum et al. |IJ and has been 
confirmed experimentally 0- This concept, introduced 
in the field of nonlinear dynamics, provides a measure 
of synchronization alternative to conventional linear ap- 
proaches. It may be useful for biological time series, in 
particular to the study of electroencephalographic (EEG) 
signals, where synchronization phenomena are expected 
to play a major role for establishing the communication 
between different regions of the brain |3(. 




12 3 4 

time (sec) 



FIG. 1: Examples of EEG signals; data are taken from a 
migraine patient subject to 9 Hz flash stimulation, and corre- 
spond to two frontal electrodes (Fl and F2). The full record 
is 40 sec long, only a 4 sec segment is shown. The signals are 
filtered in the alpha band. 

Migraine is an incapacitating disorder of neurovascular 
origin, which consists of attacks of headache, accompa- 
nied by autonomic and possibly neurological symptoms. 
It is estimated that in the USA, 5% of the general pop- 
ulation suffer at least 18 days of migraine a year, and 
more than 1% have at least one day of migraine a week 
0. In spite of a lot of research, there are still many 
unresolved issues in the pathophysiology of migraine. 
There is a tendency to believe that migraine starts with 
an underlying central nervous system disorder, which, 
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FIG. 2: Examples of spontaneous EEG; data, filtered in the 
alpha band, are taken from the same migraine patient as in 
fig. and correspond to Fl and F2 electrodes. 



when triggered by various stimuli, sets off a chain of 
neurologic and biochemical events, some of which sub- 
sequently affect the brain's vascular system. No experi- 
mental model fully explains the migraine process [5|- A 
wide range of events and conditions can alter conditions 
in the brain that bring on nerve excitation and trigger mi- 
graines. They include emotional stress, intense physical 
exertion, abrupt weather changes, flickering lights, and 
many others. The question we address here is: how does 
the response of migraine patients to such events differs 
from those of healthy persons? To address this problem, 
we investigate synchronization phenomena in Electroen- 
cephalograms (EEGs) recorded from migraine patients in 
presence of repetitive visual stimuli (steady-state visual 
evoked potentials, SVEPs 0), and study how synchro- 
nization between different brain regions varies in presence 
of external stimuli (i.e., while brain is processing exter- 
nal information). We find that migraine brains show in- 
creased alpha band phase synchronization, while healthy 
persons show a decreased one. Our results suggest that 
migraine patients have an overactive regulatory mecha- 
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nism, prone to instability, which renders them more sen- 
sitive to environmental factors. 

Our data are as follows. EEG is recorded from fifteen 
patients affected by migraine without aura |7J, in pres- 
ence of visual stimuli. During the acquisition, flash stim- 
uli are presented to the subjects repetitively at a rate of 
3-6-9-12-15-18-21-24-27 Hz. The mean age of patients is 
38.7 years (range 24-48 years). Each frequency of stimu- 
lation is delivered by a flash settled at a luminance of 0.2 
joules for at least 20 seconds; an interval of 20 seconds 
is interposed between the different trains of stimulation. 
EEG data are recorded by 18 scalp electrodes, placed ac- 
cording to the International 10-20 system, referred to CZ 
derivation. Impedance is settled below 5 Kil, EEG is 
digitally filtered off line by means of a digital filter with 
a band-pass of 0.3-30 Hz; the sampling rate is 128 Hz. 
Examples of EEG signals are shown in fig. ^ a l so sponta- 
neous EEG (i.e. in the absence of stimuli) is recorded for 
all patients, see fig. [3 AH patients are in the inter-ictal 
state, the time from the end of the last attack being at 
least 72 hours. Moreover, EEG data from fifteen healthy 
subjects (ages ranging from 22 to 45 years) are measured 
so as to have a control group. 

We recall how to detect n : m phase synchronization 
in noisy scalar signals |3| • Based on general theorems on 
analytical functions the following relation holds 
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This equation is known as the Hilbert transform and 
it is used to form, starting from a signal s(t) — Re£(i), 
the analytic signal ((t) = A{f)e l<t ' { - t \ where A(t) = 
\J s 2 (t) + s 2 (t), and s(t) = Im((t). To control the pos- 
sible synchronization of two signals Si(t), S2(t) the fol- 
lowing procedure is applied: the instantaneous phases 
4>\{t) and <f>2(t) are computed and the so called general- 
ized phase differences 
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with n, m integers, are evaluated Q. Phase synchroniza- 
tion is characterized by the appearance of peaks in the 
distribution of ^n tm and quantified by comparing the ac- 
tual distribution with a uniform one in the following way. 
The n : m synchronization index of si and S2 is defined 
ss Pn.rn — [Smax ~ Sn,m]/ S m ax, where S n ^ m is the Shan- 
non entropy of the distribution of $ n , m and S max is the 
entropy of the uniform distribution; in the case at hand 
only 1 : 1 synchronization leads to interesting results. 

Let us now turn to describe our findings. The EEG 
signals are filtered in the alpha band (8-12.5 Hz) [ljj 
and the synchronization index above described is eval- 
uated for all pairs of electrodes, for all thirty subjects 
and for all frequencies of the flash stimuli. These indexes 
are subsequently averaged over all the possible pairs of 
sensors, for each subject both in presence of stimuli and 
in spontaneous conditions. These mean values do not 
separate patients from healthy subjects; what emerges 
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FIG. 3: (a) The FDR method (see the note 0) is applied to 
select the stimuli frequencies separating patients and controls 
according to T values. The vertical axis represent the prob- 
ability that the thirty F values of patients and controls were 
drawn from the same distribution, according to the paired t- 
test. Frequencies 9-24-27 Hz are selected with false positive 
rate 0.05. (b) The probability that the thirty T values of pa- 
tients and controls were drawn from the same distribution is 
now evaluated according to the Wilcoxon rank sum test. FDR 
selects frequencies 9-24-27 Hz with false positive rate 0.05. (c) 
The FDR method is applied to select separating electrodes, 
for 9 Hz flash stimuli. The vertical axis represent the prob- 
ability that the thirty T s values of patients and controls, for 
each sensor s, were drawn from the same distribution, accord- 
ing to the paired t-test. Eleven electrodes, out of eighteen, 
are selected. The labels for electrodes correspond to Inter- 
national 10-20 system, (d) As in (c) for 24 Hz flash stimuli. 
Thirteen electrodes are selected. 



as correlated with the migraine pathology is the ratio 
7 = p{ i/ Px'x, where p{ x is the mean phase synchroniza- 
tion in presence of flash stimuli, whereas p± j is the mean 
spontaneous phase synchronization. This ratio measures 
how phase synchronization varies, in the presence of the 
stimuli, with respect to basal conditions, i.e. the neat ef- 
fect of the stimulus. Our supervised analysis (hypothesis 
testing) shows that the index T = In (7) [Tl| separates 
the class of patients and the class of controls for stimulus 
frequencies of 9,24,27 Hz. For each of the 9 flash stimuli 
frequencies u>, we apply the paired t-test to evaluate the 
probability P^ that indexes Ts were drawn from the same 
distribution (the null hypothesis); in seven cases out of 
nine this probability is less than 0.05, the standard value 
used in literature to reject the null hypothesis. How- 
ever, here we deal with multiple comparisons. To control 
the number of false positives, we use the false discov- 
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ery rate (FDR) method [Tsj. This procedure 01 selects 
the stimuli frequencies 9-24-27 Hz as separating patients 
from controls (see figEt)j with the expected fraction of 
false positive 0.05. The same frequencies (9-24-27 Hz) 
are selected by use of the standard Bonfcrroni correction 
for multiple comparisons 0] as well as by FDR if prob- 
abilities are evaluated by the non-parametric Wilcoxon 
rank sum test (see figOJ 3 )- 
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FIG. 4: In the case of 24 Hz stimuli, the mean of V a (over 
patients and over controls) is represented for all the eighteen 
electrodes. On the average, phase synchronization increases 
for patients and decreases for controls. 

A topographic analysis is also performed, in order 
to check whether this phenomenon is localized in some 
cortex region. We evaluate, for each sensor s, T s = 
ln((yo{ i) s / (pTi) s) , where (-) s means averaging only over 
the pairs where s is one the two sensors. For each 
frequency of stimuli, we apply FDR method to select, 
among the eighteen electrodes, those separating patients 
from controls according to their T s |Tfij . The results are 
depicted in Figures^ (9 Hz case) and|3Jl (24 Hz): eleven 
electrodes are recognized as separating in the case of 9 
Hz stimuli and thirteen in the case of 24 Hz; no electrode 
is found to be individually separating when 27 Hz stimuli 
are considered. Since separating electrodes from all the 
regions of the cortex (frontal, parietal, central, temporal 
and occipital) are found, it follows that the phenomenon 
here described is extended over all the cortex, not being 
localized in a limited region. Its diffuse nature suggests 
that genuine spatial synchronization UM is here involved; 
indeed, volume conduction effects [l7j would induce spa- 
tially more localized change. 

Our data show that, for patients, the mean phase 
synchronization increases in presence of visual stimuli, 
whereas it decreases in controls. For example, in the 
case of 24 Hz stimuli, and for all the sensors, the mean 
value (over subjects) of T s is shown in fig. 4: hyper phase 
synchronization is observed in patients, whereas healthy 
subjects show a reduced phase synchronization. Similar 
patterns occur for 9 and 27 Hz stimuli. In fig. 5 the his- 
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FIG. 5: Histogram of $ M for the pair T3-T5. (a) A healthy 
person without stimuli, (b) The same healthy person in pres- 
ence of 9 Hz stimuli, (c) A patient without stimuli, (d) The 
same patient in presence of 9 Hz stimuli. 



tograms of $1,1, corresponding to electrodes T3 and T5, 
are shown for a migraine patient and for a control, both 
under stimulation and spontaneously. The distribution, 
when stimuli are delivered, broadens for the healthy per- 
son while becoming more peaked for the patient. This 
behavior is further illustrated in figEl where the time 
evolution of the phase difference between two sensors is 
depicted for a migraine patient, both subject to stim- 
uli and in spontaneous conditions. In presence of flash 
phase locking, in the two signals, is observed for time seg- 
ments several seconds long; no such locking is observed 
in the spontaneous case. Phase difference curves, for a 
control, are drawn in figd It is worth stressing that this 
phenomenon is not mined if coherence is used to measure 
synchronization: considering the linear index obtained by 
integration of the coherence function (normalized ampli- 
tude of the cross spectrum of the two time series 0|) in 
the alpha band, the corresponding T and T s quantities do 
not lead to separation between patients and controls for 
any frequency of stimulation. We show that migraineurs 
are characterized by alpha band hyper-synchronization 
in presence of visual stimuli. We also show how this 
varies with the frequency of the flash, and present a to- 
pographic analysis where separating electrodes are rec- 
ognized. Whilst it is comprehensible that 9 Hz stimuli 
might cause hyper-synchronization in the alpha band (8- 
12.5 Hz), in order to figure how 24-27 Hz stimuli may act 
on alpha oscillations we observe that brain is a nonlinear 
system, and sub-harmonics of 24-27 Hz fall in the alpha 
band: stimulation in the 24-27 band may cause hyper- 
synchronization through their sub- harmonics. However 
a similar behavior is not observed for other frequencies 
with sub-harmonics in the alpha band, like 18 Hz: fur- 
ther investigation is needed to clarify this aspect of the 
phenomenon. It will be also interesting to investigate the 
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FIG. 6: Time evolution of A$ (i.e., the phase difference 0i — 
4>2, without wrapping in the interval [— tt, tt]) for the pair Fl- 
F2 in a migraine patient (the same patient as in fig. |5Jl. Top: 
in presence of 9 Hz stimuli. Bottom: without stimulation. 

response of migraine patients with aura. Our results are 
consistent with current theories about the role of subcor- 
tical structures in migraine. Since brainstem is active in 
migraine [l9j , it has been proposed, as a unifying concept 
of migraine, that brainstem regions concerned with neu- 
ral mechanism of synchrony are dysfunctional |20| . Cor- 



tex, in migraine brains, is thus misled by a dysfunctional 
gating system; normal light is unpleasant, normal sound 
uncomfortable and, probably, normal pulsing of vessels 
felt as pain. On the mathematical side, our results con- 
firm the usefulness of the analytic signal technique to 
study physiological time series. 
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FIG. 7: Time evolution of AO for pair F1-F2 for a control. 
Top: in presence of 9 Hz stimuli. Bottom: without stimula- 
tion. 
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